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Abstract. Motivated by some numerical observations on molecular dynamics simula- 
tions, we analyze metastable trajectories in a very simple setting, namely paths generated 
by a one-dimensional overdamped Langevin equation for a double well potential. More 
precisely, we are interested in so-called reactive paths, namely trajectories which leave 
definitely one well and reach the other one. The aim of this paper is to precisely analyze 
the distribution of the lengths of reactive paths in the limit of small temperature, and to 
compare the theoretical results to numerical results obtained by a Monte Carlo method, 
namely the multi- level splitting approach [6]. 



1.1. Motivation and presentation of reactive paths. A prototypical example of a 
dynamics which is used to describe the evolution of a molecular system is the so-called 
overdamped Langevin dynamics: 



where x\ 6 R d denotes the position of the particles (think of the nuclei of a molecule) , 
V : R d -)• R is the given potential function modeling the interaction between the parti- 
cles, (Bt) t>0 is a standard Brownian motion on R d and e is a (small) positive parameter 
proportional to temperature. The potential V is assumed to be smooth and to grow suf- 
ficiently fast to infinity at infinity so that the stochastic differential equation ([TJ admits a 
unique strong solution. One common feature of many molecular dynamics simulations is 

that the dynamics ([I]) is metastable: the stochastic process (^x[ £ ^ spends a lot of time 

in some region before hopping to another region. These hopping events are exactly those 
of interest, since they are associated to large changes of conformations of the molecular 
system, which can be seen at the macroscopic level. 

In the following, we focus on the limit of small temperature (namely e goes to zero). 
In this case, the Freidlin-Wentzell theory |11| is very useful to understand these hopping 
events. Specifically, it turns out that the metastable states are neighborhoods of the local 
minima of the potential V, and that the time it takes to leave a metastable state to reach 
another one is of the order of 



Here, 5V is the height of the barrier to be overcome (namely the difference in energy 
between the saddle point and the initial local minimum), and C is a constant depending 
on the eigenvalues of the Hessian of the potential at the minimum and at the saddle point 
(see Equation ([3| below for a precise formula in the one-dimensional case). This is the 
so-called Eyring-Kramers (or Arrhenius) law, and we refer for example to U [2j [15] for 
more precise results. 

Actually, the most interesting part of a transition path between two metastable states 
is the final part, namely the piece of the trajectory which definitely leaves the initial 
metastable state and then goes to the next metastable region: this is the so-called reactive 
trajectory (or reactive path) 1 12^ 110]. In particular, reactive paths give important informa- 
tion on the transition states between the two metastable states. One numerical challenge 
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in molecular dynamics is thus to be able to efficiently sample these reactive paths. Notice 
that from the Eyring-Kramers law pi), a naive Monte Carlo method (generating trajecto- 
ries according to ([I]) and waiting for a transition event) cannot provide efficiently a large 
sample of reactive paths, hence the need for dedicated algorithms. 

In [6], we proposed a numerical method based on an adaptive multilevel splitting al- 
gorithm to sample reactive trajectories. One interesting observation we made is that the 
lengths of these reactive paths seem to behave very differently from Q, see Figure [l] below. 
It seems that, in the limit of small e, the distribution of these lengths is a fixed distribu- 
tion shifted by an additive factor — log e. The aim of this work is to use analytical tools 
to precisely analyze this distribution in the asymptotic regime e goes to zero, and to give 
a proof of this numerical observation. 

1.2. The one-dimensional setting and our main results. In the following, we con- 
sider a one-dimensional case (d = 1), and we assume (for simplicity) that the potential 
V admits exactly two local minima (V is a double- well potential). More precisely, let us 
denote x* < y* the two local minima of V and z* £ (x* , y*) the point where V reaches its 
local maximum in between. As explained above, we are interested in trajectories solution 
to ([I]) from x* to y* , and more precisely in the end of the path from x* to y* (the reactive 
paths). In order to precisely define these reactive paths, let us introduce the first hitting 
time of a ball centered at y* with (small) radius 5 y > 0, starting from x*: 

Ty* = inf it > : \x\ e) - y*\ < 6 y } with X ( Q £) = x* . 

In this setting, formula Q writes (notice that V"(x*) > and V"(z*) < 0): 

( 3 ) E ( T y : ) ~ n ft r iii 2 Lnu W 6Xp ~ V ^/ e ) • 

V y / ^V"(x*)\V"(z*)\ 

The cf-dimensional version of this result is established in [J]. Let us also introduce the 
last exit time from the ball centered at x* with (small) radius 5 X > before the time Ty* 

(again starting from X^ = x*): 

Sf* = sup {i < : \X<f ] - x*\ < 5 x y 

The question we would like to address is: how long is a reactive path, that is the time 
T*: - Sf* as e -> ? 

This question was partially addressed in jTT] where the ball centered around y* is re- 
placed by the complementary of the domain of attraction of x* for the deterministic dy- 
namical system corresponding to ([T]) with e = 0. Several papers are dedicated to the more 
subtle situation where points on the boundary of this domain are not attracted to x*. In 
our simple framework, such a domain is given by (— oo,z*) (see |15| for such a study). In 
19], Day is interested in the law of the exit time from a domain containing an unstable 
equilibrium when the diffusion starts on the stable manifold. Thus, even if the laws of the 
exit times considered in these papers are related to the distribution of the lengths of reac- 
tive paths we deal with in the present work, these are different quantities, with different 
asymptotic behaviors. 

In order to specify our purpose, let us now make our assumptions on the potential V 
more precise. 

Assumption 1.1. The potential V is smooth, has exactly two local minima x* < and 
y* > and a local maximum z* = 0. Moreover, V' is positive on (a;*,0) and negative on 
(0,y*) and the local maximum at is assumed to be non-degenerate: 

(4) V(0) = 0, V'(0) = 0, and V"(0) = -a < 0. 
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Notice that the potential V is close to x i— )■ —ax 2 /2 for values of x around 0. More 
precisely, it is easy to show that there exist K > and 5 > such that, for all \x\ < 5, 

(5) -ax-ifx 2 < V'(x) < -ax + Kx 2 and — < V(x) < 1 — . 

Example 1.2. An example of a potential which satisfies the assumption 1.1 is 

4 2 

(6) V: x^—-—. 

In this case, —1 and +1 are £/ie two (global) minima. This is a double well potential with 
a local maximum at x = which is non degenerate, with a = 1. 

Let us denote A = x* + 5 X G (x*,0), B = y* - 5 y £ (0,y*) and x G (A, 0). We are 
interested in the behavior of 

T x ->b = inf |i > : = i?| conditionally to the event |-^o^ = x ' < Ta\ 

when e goes to zero. At the end of the day, the aim is to let x go to A. As mentioned 
above, simulations in [6] suggest that, if the local maximum is non degenerated, then the 
law of this length looks like a fixed law shifted as e goes to 0. Figure [T] presents the density 
of the reactive path T x ^b f° r several values of e, when V(x) = x 4 /4 — x 2 /2, A = —0.9, 
B = 0.9, and x = —0.89. In (6j E], it is suggested that the asymptotic shape of these laws 
is an Inverse Gaussian distribution. In fact, it is not the case: it turns out to be a Gumbel 
distribution. 

Definition 1.3 (Standard Gumbel distribution). The standard Gumbel distribution is de- 
fined by its density function 

f(x) = exp (—x — e~ x ). 

Its Laplace transform is given by 

sG\ 



E(e- sG ) 



r(i + s) ifs>-i, 

+oo otherwise, 

where T(z) = J °° t z ~ 1 e~ t dt is the Euler's Gamma function. 

The main result of the paper is the following convergence in distribution. 

Theorem 1.4. Under Assumption \1.1\ for any A G (x*,0), B G (0,y*) ; and x G (A, 0) 
we have, conditionally to the event |Aq = x, Tb < Tj^, 

T X ^ B + ~ loge -(log(|x|5) + F(x) + F{B) - loga + G) 

a e->o a 

where G is a standard Gumbel random variable and 

Fis) =J°{m +l t ]dt 

for any s G (x*, y*). 



Notice that by the integral defining the function F is well defined. We slightly abuse 
notation and denote Ta^b the limit of T x ^b when x goes to A. We then have 

Ta^b + - loge -{\og(\A\B) + F(A) + F(B) - loga + G). 

a e-s>o a 

Example 1.5. Let us come back to our previous example where the potential V is defined 
as 

A 2 
X X 

V : x h-> . 

4 2 
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Figure 1. Left : Density of the length T x ^b f° r different values of e (from 
left to right, e = 1,0.5,0.2,0.1,0.05,0.02,0.01) when V(x) = x 4 /A-x 2 /2, 
A = —0.9, B = 0.9, and x = —0.89. Right : Empirically centered versions 
of these densities. 



In this case, a = 1 and if we choose A = —0.9, B = 0.9, and x = —0.89, we get 

T_ .89-,o.9 + loge log(0.89 x 0.9) - J log(l - 0.89 2 ) - \ log(l - 0.9 2 ) + G. 
This is illustrated on the left hand side of Figure [i] and on Figure [3] below. 

The paper is organized as follows. Section [2] recalls classical tools that are used in the 
proofs. Section [3] provides a key estimate for the (repulsive) Ornstein-Uhlenbeck process. 
The proof of Theorem 1.4 is given in Section [4] Finally, Section [5] is devoted to particular 
potentials that are degenerated at the origin {i.e. V"(0) = 0) or singular {e.g. V{x) = 
-|s|). 

2. Classical tools 

2.1. Laplace transform of the exit time. Let us first recall how one can link the 
Laplace transform of the exit time of an interval to the infinitesimal generator A £ of the 
diffusion process ([T]) where 

A £ f{x) = ef"{x)-V'{x)f l {x). 
Fix a < x < b and denote by H^l the first exit time from (a, b), starting from x: 

= inf jt > : x[ £) £ {a, 6)} = A T b (e) where = inf jt > : x[ s) = c}. 

In the sequel, we may drop the superscript e and the indices a and b and simply denote H 
for H { a % 

Notice that = b\ = {T^ < T a }. For any s G [0, +oo) and x G (a, b) let us define 

(7) F £ (s,x) := E x (e- sH |xg } = ft) and F £ {s) = lim F £ {s,x). 
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Let us also introduce the function u, solution of 



(8) 



A e u s (x) = su s (x), x £ (a, b), 
u s (a) = 0, u s (b) = 1. 



Ito's formula ensures that {u s {X^)e st ) t>0 is a martingale and then 



E x (u s (X$)e- sH ) =E x (e~ sH Ir „ = A 



Consequently, 

(9) F £ (s,x) = ^\. 

This formula will play a crucial role in the following. 

Remark 2.1 (The exit distribution). When s = 0, Equation ([8j is easy to solve: for any 
x G (a, 6), 

r* e v{,)/ e ds 

J a 



u Q (x)=F x (T b <T a ) 



r b e v(s)/e ds 

J a 



2.2. The /i-transform of Doob. The process (jf ) i>0 solution of the stochastic differ- 
ential equation ([T]) conditionally to the event {Tt, < T a } is still a Markov process. Moreover, 
it can be seen as the solution of a modified stochastic differential equation with a drift that 
depends on the exit probabilities for the process. This is the so-called h-transform. 

Proposition 2.2. Conditionally to the event {Tf, < T a }, the process X^ E ' is a diffusion 
process and it is the solution of 

(10) dX ( t £) = VTedB t + ( -V'(X^) + 2s f ^Ql {Tb>t} ) dt 



where, for any x £ (a, b), 

f x e v ^/ £ ds 

J a 



V h £ {X\ 
r x p V{s)/e n 

h e (x) 



r e V(s)/e ds 
J a 

See [8] for the proof of this assertion via Girsanov's theorem. Similarly, one could write 
the equation satisfied by a diffusion process conditioned to reach a given point at a given 
time (see [16J for instance). 

Remark 2.3. The additional drift is singular at point a and is equivalent to 2e(x — a) -1 . 
This ensures that Y cannot hit a as far as t < T a (see the Feller condition in [17\). 

Let us associate to a potential V the modified drift induced by the /i-transform on the 
interval (a, b): 

h! (x) e V(*)/£ 

(11) b v (x) = -V'(x) + 2e-4^ = -V'{x) + 2e- 



Lemma 2.4. Let us assume that x*<a<0<b<y* and that V satisfies the assump- 
tion \l.l\ Then, for any x £ (a, b), 

b v (x) — > \V'(x)\. 

Proof. Since V is increasing on (a, 0) then, for any x £ (a, 0), 

p V(x)/e 

i v ^/ £ ds ~ e—— and b v (x) ~ V (x) = \V (x)\ . 

e^0 V'(x) e-5-0 1 1 
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In other words, the /i-transform turns the negative drift —V'{x) to its opposite. Moreover, 
it is obvious that, for any x > 0, h e (x) goes to 1 as e — > and h' e {x)/h e {x) goes to 
exponentially fast: in this case, by(x) —> —V'{x) = \V'{x)\. Finally, one can notice that 



since V(s) ~ V" (0)s 2 /2 when s goes to zero. □ 

The /i-transform and the previous Lemma will be two major ingredients for the argu- 
ments below. 

In the former proof, and in the following, we constantly use the Laplace's method to get 
equivalents of integrals in the limit e goes to 0. Let us recall these classical results: 

Lemma 2.5. Let [a, 6) be some interval of K (with possibly b = oo), tp : [a, b) — > M a 
function continuous at point a such that ip{a) ^ and tp : [a, b) — > K a function of class C 2 



such that <p' < on (a,b). Let us denote /(e) = exp((p(x)/£) ip{x) dx. Then, we h 
• If<p'(a) = and ip"(a) < 0, 



ave: 



lf<f/(a)<0, 



/ TT£ 

f{£) sZo^^W\ eXMa)/£) ^ a) - 



f( £ "> ~n I if \ I ex P(¥ : '( Q )/ £ )V ; (a)- 
e-s>o \tp (a)\ 



3. Main example: the repulsive Ornstein-Uhlenbeck process 

In this section we deal with the simplest example of a potential that is smooth and 
strictly concave at the origin. We assume here that V(x) = —ax 2 /2 on the set [—6, b] with 
b, a > and then investigate the behavior of the process: 



(12) dYf' = V2edB t + aY^> dt. 

In the sequel we denote 

T^' a,x = inf {t > : Y t (£ ' a) = &} with Y^> a) =x£ (-6, b). 

For the sake of simplicity, we first deal with the case a = 1 and then we will get the 
general result thanks to a straightforward scaling. The strategy is to express the Laplace 
transform of this exit time in terms of special functions and then to derive its asymptotic 
form as e goes to 0. In the sequel, Tj, stands for T^ £ ' l,x \ 

Proposition 3.1. Let x £ (—6,6). For any s > —1, we have 

'r(l + s ) e -«(-log e +log6+log| a: |) i j x£ (_5 i0 ), 

2 s /2 /l + s N 



(13) i a: (e- s ^|r 6 <r_ 6 



E-S>0 



A \ 2 / 



e -s(log6-logx) ifxe{0,b). 

One can also notice that lim e _s.o (e~ sTb | Tj, < T_b) = oo if s < — 1. 
Proof. The Laplace transform of the exit time is linked by ^ to the solution u s of 

{£u'g(x) + xu' s (x) = su s (x), x € (—6, 6), 
u s (-b) = 0, 
u s (b) = 1. 
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Let us define b e = bj^fe and the function v s on (—b e ,b e ) by v s (y) = u s {yyfe). Then v s is 
the solution of 



(15) 



v'siv) + yv' s {y) = sv s (y), y G {-b £ , b E ), 
v,(-b e ) = 0, 

[v.s(b £ ) = 1. 



As it is recalled in Section 2.1 (see (|9j)), one has 

E,(e-^|T 6 < T_ 6 ) = *W = V -fl^. 

u (x) v (x/^e) 

One can express the function v s in terms of some special functions. Let u > and define 
the parabolic cylinder function D- v as 



DJx) 



1 f'°° 

r O) Vo 



The so-called Whittaker function -D_„ is solution of 

Z^(*)- ^ + I /-^D_,(x) = 0. 

See pQ ch.19] or j3j p. 639] for further details. Define the function ip u by 

Vu (x) = e- x2 / A D^(x). 
One can check with a straightforward computation that 
(16) <p"(x) + xip' v {x) = (y - l)<p v (x). 

In the sequel, s and v are linked by the relation 

v = s + 1 > 0. 



Notice that £ i— )■ ip v {—x) is also solution of (16) (and ip u and y„ are linearly inde- 
pendent). Then, the solution of (15) is a linear combination of (p u and tp v satisfying the 
boundary conditions. The function v s is given by 



(17) 



v s (x) 



Let us study the asymptotic behavior of (p u (b) and (p u (—b) as b — > +oo. The Laplace's 
method ensures that 

+2/ .. r(i/) 



t u-l e -t-/2 e -U dt 







As a consequence, 



Moreover, 



b-^+oo b h 



<Pu(-b) 



t ,-l e -(t-bf/2 dt 



'2ir. 



M-l 



In particular, one obtains that 

M-b) 2 ~ Mb) 2 



b-t+oo 



M-bf 



2tt 



b-^+oo r(z/) 2 



&2(>-l) 
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Moreover, we get, for any 7 £ (0, 1), that 



<p v (-b)(p u (yb) - (-76) 



e -7 2 b 2 /2 



b-s>+oo T(u) (76) 1 



V2tt 



(76) 



-6 2 /2 



6->+oo r(z/) 7^6 

As a conclusion 

M-b)M7b)-Mb)M-lb) T{ v )e-^ b2 ' 2 I» e^ 2 ' 2 



(18) 



¥v{-b) 2 - tp v (b)< 



u-X ■ 



One can then deduce the asymptotic behavior of v s solution of Equation ( 15 ) at the point 
x/y/e (with x < 0) replacing in Equation (18) b by b £ = b/y/e and 7 by —x/b with 
7 6 (0, 1). Since v = s + 1, this leads to 



V s (x/y/e) 



r(i + s ) 



-z 2 /(2e) 



and 



^0 (-^/^) s+1 (Vv^) 
r(l + a) 



v s (x/y/e) ^ 



Tl + s) 



e 

x\b 



This is the expression of the Laplace transform in Equation (13) when x S (—6,0). The 
two other cases are easier to deal with. If x = 0, one has (since vq(0) = 1/2) 

v s (o) 2^(0) /21 r +0 ° t s e -ty2 dt = ^ r (i + Ai 



Vq(0) e^0 ip u {-b £ ) e^0 \ irb s £ J 
At last, if x = 76 with 7 £ (0, 1) then 

u s (7& £ ) ip v {-ib e ) 
^0 (tM £ -^° fiy(-be) 



□ 



Remark 3.2. T/ie parabolic cylinder functions D- u also appear in [5j, Section 2, where 
the author studies the first exit time from a square root boundary for the Brownian motion. 

Proposition |3.1| yields the following convergence in distribution. 

Theorem 3.3. Let a > and x 6 (—6, 6). Conditionally to the event ^T^ e ' a,x ^ < T^ a ' x ^\, 
we have 

— loge + log(|x|6) + G — logo if x G (—6, 0), 

— log y/e + log b + G — log y^ct if x = 0, 
log 6 — logs if x £ (0, b), 



(19) 



-,(e,a,a;) C 1 



e->0 a 



where the law of G is the standard Gumbel distribution and G is a random variable with 
Laplace transform given by 







E e 



-sG 



2 s / 2 fl + s 



'7T 
+ OO 



if a > -1, 
otherwise. 



Proof. The case a = 1 is a straightforward consequence of Proposition 3.1 Moreover, for 
any positive constants r and a one has, for any t > 

aY^'^ = aY (£ ' a) + aV2eB Tt + <ra / ds 

Jo 

= aY (e ' a) + Wo^V^Bt + «t / aY T ^ du. 

Jo 
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FIGURE 2. Mean length of the reactive path for the repulsive Ornstein- 
Uhlenbeck process dY t (£) = y/2edB t + Y t {e) dt, with Y (e) = -0.89, on the 
set [—0.9, 0.9] as a function of log e (see Equation ( 19 )). The 95% confidence 
intervals are of the size of the points. The function log e i— > — log e + log(| — 
0.89| x 0.9) + 7 is drawn in dotted line. These results have been obtained 
with the algorithm described in [6]. 



This ensures that if a = \fa and r = 1/a, then the process (cly 



is solution of 



Equation (12) with a = 1 and the initial condition uYq . In particular 



£^(rp{s,a,x) I rp(e,a,x) rp(e,a,x)\ £ / —lrp(s,l,x/y/a) i rp(e,l,x/^/a) rp{e ,1 ,x / y/a) 

\ b \b —b J \ b/y/a ' b/y/ot -b/y/a 

The result for a 7^ 1 is then a straightforward consequence of the result for a = 1. 



□ 



Notice that the formulas (19) admit a limit when x goes to —b. Before coming back 



to the general case, let us conclude this section with a few remarks about the case of the 
Ornstein-Uhlenbeck process. 



Remark 3.4. Let us discuss the asymptotic behavior (19) of the length of the reactive path 
when x £ (—6,0) and e goes to 0, taking for simplicity a = 1. The time \og{b/\/e) is the 
time needed by the deterministic process y( 0,1 ) to go from \fe to b since y^ ' 1 ) = e i y/e. The 
Freidlin-Wentzell theory tells us that the first part of the reactive path (from x to —\fe) has 
a similar length \og(\x\ / \fe) . Finally, the Gumbel variable G accounts for the (asymptotic) 
random time needed by Y^' 1 ) to go from —\fe to \fe. 

Remark 3.5. It is easy to check from the proof that the results of Proposition [XT] are still 



valid if b = b e and x = x e depend on e as long as b e j \fe and x £ j \fe go to infinity when 



£ goes to zero. For example, if b £ > is such that lim £ _>.o b £ /\/e = oo and x £ 6 (— b £ ,0) 
is such that lim £ ^.QX £ /\/e = —oo, then r r^ e ' a,Xe ' > ~^ i (— log e + log 
This remark will be useful in SectionV 



+ G-\oga). 



Remark 3.6. Figure^ illustrates Th 



eorem 



cess dY, 



2edB t + Y t {£) dt, with Y (e) 



3.3 for the repulsive Ornstein-Uhlenbeck pro- 
= -0.89, on the set [-0.9,0.9]. Denoting 



T_o.89->o.9 the length of the reactive path from —0.89 to 0.9, then Equation (19) ensures 
that E[T_o.89->0.9] i- s equivalent to — logs + log(| — 0.89| x 0.9) + 7, when e goes to zero 
(7 stands here for the Euler's constant). Figure^ compares this theoretical result with the 
empirical means obtained thanks to the algorithm described in [6] for e ranging from 0.01 
to 1. 
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4. The general (strictly convex) case 

Let us now come back to the general strictly convex case described in Section [T] We 
recall the notation. The potential V has exactly two local minima x* < and y* > and 
a local maximum z* = 0. Moreover, V is positive on (x*,0) and negative on (0, y*) and 

V(0) = 0, V'(0) = 0, and V"(0) = -a < 0. 

Let us consider A G (x*,0), B G (0, y*) and x G (^4,0). We are interested in the behavior 
of 

T x ^b = hif (t > : -X"^ = i?| conditionally to the event |-Xq £ ^ = x > < Ta\ 
when e goes to zero. 

According to the Markov property, and considering the initial point x G (^4,0), the 
strategy is to decompose the reactive path from x to B into three independent pieces: 

(20) H = T x ^.^ Ce + T_ C£ ^ 6e + T b£ ^ B 

on the event {Tg < Ta} where < c £ < b £ < \x\ A B will be chosen in the sequel. More 
precisely, we will choose 

a 2 1 

b £ = e p and c £ = e 1 with - < (3 < 7 < — . 

5 2 



The first and third times in (20) are essentially deterministic, as specified by the following 
result. 

Proposition 4.1. If < /3,7 < 1/2, then, conditionally to the event {Tb < Ta}, 

p p 

Tb e ^B — h e ~^B > and T x ^_ Ce — t_ Ce -kr > 0, 

e— >0 e— »0 

where tb £ ->B is the time for the unnoised process to reach B from b £ G (0, B): 

' KB = - / T77TT ds, 



V'(s) 

and t-c^x is the time for the unnoised process to reach x from —c £ G (x, 0): 

1 y r (s) ds ' 



This is proved in Sections 4.1 and 4.2 In Section |4J3| we compare the second time in (20 ) 
to the reactive time of an Ornstein-Uhlenbeck process. 

4.1. Going down is easy. The easiest part is to study the third time interval T{, £ ^b- 
Our goal here is to prove that, starting at b £ , the process is close to the deterministic 
path (xt) t>0 solution of the ordinary differential equation 

(± t = -V'(x t ) t>0, 
\x = b £ . 

In this aim, we need to state a few intermediate results. First, it is readily seen that, 
starting at b £ , the probability for the process (Xj f )t>o to hit before B goes to at an 
exponential rate when e goes to 0. Indeed, since b £ = e 13 with f3 < 1/2, we have 

f h B e v ^ds I 2ae v(h ,, 

bA ° B) f B eV(^ds ^0 y n\V'(b e )\ 

In the following, we will denote Sl £ the event on which this does not occur, so that P(O e ) 
goes to 1 when e goes to 0. 

Of course, this will also be true for the event Q, x which is defined as: the process starts 
at a fixed point x G (0, B) (independent of e) and does not hit before B. Again, P(JXe) 
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goes to 1 when e goes to 0. Then, starting at x G (0,B), our aim is to compare the 
deterministic path (xt) t>0 solution of the ordinary differential equation 



(21) 

and the random process 



x t = -V'(x t ) t>0, 

X = X, 



'dX ( t s) = -V'{X[ 6) ) dt + V2edB t t > 
= x. 



For this, let us introduce c G (B, y*) such that c — B < B — x, the deterministic time t c 
t x ^c = m f {t > : Xt = c} and the stochastic time T c = T x ^ c = inf It > : X^ = ( 



Lemma 4.2. Define K := sup s6 [ 0)C ] |V"(s)|, then 

P(n x n< sup |xj e) -x s \ > 7] \ ) < 2exp ( - 

\ 0<s<t c AT c I \ 



^y2g ^Ktc 

4et c 

Proof. Let us assume that we work on the event Q x . For any t < t c A T c , 



X ( t £) -x t = - f\v'(X^) - V'(xs)) ds + VteBt. 
Jo 



The Gronwall Lemma ensures that 



sup \X( £) - x s \ < V2ee K{tcATc) sup \B S \ < V2ee Ktc sup \B S \. 

0<s<t c AT c 0<s<t c AT c 0<s<t c 

Finally, the reflection principle for the Brownian motion ensures that sup 0<s<t B s has the 
law of \Bt\. As a consequence, for any r, t > 0, 

P ( sup |£ s | > r ] < 2P ( sup B s > r ) = 2F(\B t \ > r) = A¥(B t > r) < 2e~ r ' l{2t) . 
\0<s<t ) \0<s<t ) 

This concludes the proof. □ 

The first consequence of this result is that the stochastic time T X ^ B required by the 
random process to go from x G (0, B) to B converges to the deterministic time t x ^ B as 

£ -> 0. 

Corollary 4.3. Let < x < B, then 

P 

T X ->B -> t x ->B- 

Proof. Since we will apply the result of the previous lemma, we still work on the event 
Q x . Let us denote 77 a real number such that < r] < c — B. Then, on the event 

\xi £) - x s \ < r/}, the random time T X ^ B belongs to the deterministic 
interval [t x -^B-ri,tx-*-B+Ti]- In other words, 

f B ds f B+r > ds 

Jb-.v 7 ^) ~ ~ tB -^ B ^ T ^ B ~ ^ B ^ tB ^ ~ ~ J B vw' 

As a consequence, for any r\ G (0, c — B), 

\T X ^ B ~ *x-kb| <V X sup 



se[B- v ,B+ v } \ v '( s )\ 
Finally, for any i] G (0, c — B), 

F(n x r)\\T x ^ B -t x ^B\>VX sup —j— \ ) < 2expf- 7? 



s£[B—r),B+r)] \V'{S)\ J J ~ ' V 4e *. 
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where t c = t x ^ c . This concludes the proof of the corollary. □ 

Our next goal is to prove that this result still holds if the starting point, namely b £ = e^, 
goes to sufficiently slowly as e — > 0, that means if /3 < 1/2. Let us fix D 6 (b £ ,B) (for 
sufficiently small e) such that 

sup \V"{s)\ < — . 

s<=[0,D] Z P 

This is always possible since (3 < 1/2, V"(0) = a, and V is assumed to be smooth. Then, 
as previously, we fix c S (D, B) such that c — D < D — b £ , and 

K:= sup \V"{s)\<~ 

s£[0,c] Z P 

Corollary 4.4. If < b e = e 13 with (3 < 1/2, then 

\Tb e -+D — tb e ^D\ — -> 0. 

Proof. Here, we work on the event Q £ , which is not a problem since, as mentioned above, 
P(f2 e ) goes to 1 as e goes to zero. The first part of the proof is similar to the ones of 
Lemma |4. 2 1 and Corollary 4.3 For any rj G (0, c — D), 

F I ft £ n | \T b ^ D - t bE _ D \ > V x sup — ^- i | < 2 exp ( ' 



sE l D - v>D+v ] \V'(s)\ J J \ 4et 

where t c = i& e _» c . Moreover, since V'(s) ~ s -»o —cts, we have 



As a consequence, 



r ° ds 1 
t b ^ c = -l — = --lo g b £ + £ (l) 



-2Kt 

e A ^ tc as oc 



~ T7T-, > +00. 



Aet c £-s>o —4/3 log e e^o 
This proves the convergence in probability of T be ^D as £ goes to 0. □ 

Finally, according to the Markov property, we can summary the previous results by 
decomposing the path from b £ to B into two independent pieces: 



Using Corollaries |4.3| and 4.4, we immediately get the following proposition. 



B ds 



Proposition 4.5. If < b £ = e 13 with (3 < 1/2, then 

IP f 
\Tb E ^B ~ h e -+ B \ where tb E ^B = - J 

Remark 4.6. IfV is given by (|6j ; one can compute the expression of the solution (xt) t>Q 
of (21). Let us define the function ^ on (0, 1) by 

ty(x) = log 

Notice that 



r,u \ 1 1/2 1/2 1 

* (x) = i " ~TT7 = "tTTTT- 

x x — 1 x + 1 V [xj 

As a consequence, the derivative ofti-t *&(xt) is equal to 1 and 

x t = ^- 1 (^(x)+t). 
Moreover the elapsed time from x £ (0, B) to B £ (0, 1) is given by 

t x ^ B = *(B) - *(x) = - log(x) + *(B) + X - log(l - x 2 ) 
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As was just proved, this result still holds when x = b £ as long as < /3 < 1/2. 
4.2. The climbing period. 

Proposition 4.7. If c £ = e 7 with 7 < 1/2, then conditionally to the event {T_ Ce < T A }, 
and for x £ (A, 0), 



Tr,->-o, ~ t-,^ T \ > where i_ r 



ds 



Proof. One has to consider the /i-transformed process and use the fact that the new drift 



converges to V'(s) uniformly on [A + 5, —c £ ] with small 5 as e goes to 0, see Lemma 2.4 



above. □ 

4.3. Central behavior. Let us finally study the behavior of T_ c ^ bs conditionally to the 
event {T be <T A }. 

The sketch of proof is as follows: 

(1) Prove that one may assume that the process does not go below — b e ; 

(2) Rescale space to change (—b e ,b e ) to (—1, 1); 

(3) Consider the /i-transform process to get the evolution of the process conditioned 
on {T x < T_!}; 

(4) Introduce the /i-transformed repulsive Ornstein-Uhlenbeck process; 

(5) Compare the drifts; 

(6) Use Theorem 4.3 in [8j; 

(7) Conclude. 

Step 1. The first step is to notice that it is equivalent to look at T_ Ce _j.{, s conditionally to 
{T be < T_ be } or conditionally to {T be < Ta}- 

Lemma 4.8. If < f3 < 7 < 1/2, there exists a constant C > such that, for any s > 

1 - an-' < E 7 I7i - <r ^ < 1 + *■>->. 

- E_ c ,(e-« H -»..».|T t , <T_ t .) - 

Proof. By continuity 

{T be < T A } = {T he < T_ 6 J U {T_ be < T bE < T A } 

where the two sets on the right hand side are disjoints. Moreover, the strong Markov 
property ensures that, for any s > 0, 



0<E_ C£ (e~ sH ^t {T _ b£<Tbe<TA} 

<m_ bs (e- sH ^t {Tbe<TA} ) 

< P-6 e (T_ Ce < T A )M_ Ce (e- sH ^t {Tbe<TA} 
As a consequence, for any s > 0, 



1 < v± 1 be<TA]j N < 1 + P- 6e (T_ C6 < T A ) < 



Making s = in this equation leads to 

TA < — ^ 

xn e <T A 



6e (T_ Ce <r A )< ^%^l <i. 
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Consequently 



l-P_ 6e (T_ C£ <T A ) < 



E_ Ce (e- sH ^\T be <T A ) 



-sH. 



\T bs < T_ be ) 



< l + P_ fe£ (T_ Ce <T A ). 



To conclude, one just has to remark that, since V( — b e ) < V(—c £ ), 

P_ 6 (T_ c < T A ) = ~ £ 7-/yv(-M-n-c £ ))/ £ < £ 7-f 



□ 



Corollary 4.9. If < (3 < 7 < 1/2, i/jen, starting from —c £ , 

C(H AA \T be < T A ) ~ q C(H Abe \T be < T_ 6 J. 

Step 2. Let us define r) £ = e/b\ and the process Y by Yi = x[ £ ^ /b £ (dropping for simplicity 

the explicit dependence on e in the notation for Y). Obviously, if is equal to — c £ , 
then Y is solution of 

l2r) £ dB t di, 



Y = 



-c £ /b £ . 



In terms of Y, we are interested in the hitting time of 1 conditionally to the event 
{Ti < 

Step 3. Thanks to the /i-transform of Doob, one can see Y, conditionally to the event 
{Ti < T_i}, as a diffusion process. Define, for any y G (— 1, 1), 



,V(b e y)/e 



J^ ie V{bss)/e ds h £ (y) J\ e nbes)/e ds - 

Conditionally to {Ti < T_i}, the process Y is solution of 

V'(b £ Y t ) h' £ (Y t ) \ 

-^ £ — + 2 ^hM ) 1 ^})dt, 



dY t = ^2r] £ dB t + 
Y = -c £ /b £ . 

Step 4. Similarly, the repulsive Ornstein-Uhlenbeck process (Zt) t>0 



I dZ t = y/2rf £ dB t + aZ t dt, 

\Z = -Ce/be, 
evolves, conditionally to the event {Ti < T_i}, as 

9'e(Zt 



where 



dZ t = ^2 V£ dB t +[aZ t + 2r] £ ^}-^t {t<Tl} 
V 9e{Zt) x ~ 

Z G = -Ce/be, 



dt, 



1 \ Ae—V^) cfa ^( y ) 
9e(y) = — T77— — and 



-ay 2 /(2 Ve ) 



f^ e -^ 2 /(2%) ds 



9e(y) J\e-^ 2 /(^e) ds' 



Step 5. Let us now notice that the drifts of the stochastic differential equations that drive 
Y and Z conditionally to the event {Ti < T_i} are close. 
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Lemma 4.10. Under Assumption [7T7J if 4/9 < j3 < 1/2, th 



en 



(22) — x sup 

Vs j/6(-l,-l] 



, 2 ^(y) 



«y + 27?, 



^(y) 
9e{y) 



£->0 



-> 0. 



Proof. Thanks to Assumption as soon as b £ < 5, we have, for any y £ [—1,1], 

V'(b e y) 



ay 



< Kb £ y 2 < Kb £ 



so that 



— x sup 

Vs y e(-i-i] 



V'(b e y) 



ay 



£->0 



-> 



as soon as /3 > 1/3. It remains to prove that sup ye (_ 1 _ 1 ] |A e (y)| goes to zero when e goes 
to zero, where 



A e (y) := 



h'M 9' £ (y) 



,V(b e y)/e 



-ay 2 /(2rj e ) 



h £ (y) 9e(y) ^eW*')/' ds f^e-^V^ds' 

We propose to do this in two steps: first for y S [—1 + e K , 1], then for y € (— 1, —1 + e K ] 
where k = (3/2. 

(1) y € [—1 + s K , 1]: Thanks to assumption |1,1[ we have for all s S [—1, 1] 



V(b £ s) = -ab 2 s 2 /2 + 6 e (s)b%g 



3„3 



with 



sup |# e (s)| < - sup V {3) (x) < C, 
se[-i,i] b xe[-i, E ,y 

where C is a constant independent of e. As a consequence, since b £ = and 

n £ = e 1 - 213 , 

e V(b e s)/e = e -c tS y(2 VE ) e 8 E (s)s3 £ 3P-^ 

Now, we can write 



e 8 ^ 3 ^ 1 = l + 5 £ (s)0 £ (s)s i e 



3^3/3-1 



with 



sup |<5 £ (s)| < e^ 3 "" 1 < C, 



«e[-l,i] 

where C is a constant independent of e. For the sake of simplicity, we denote 9 £ (s) 
for S £ (s)8 £ (s). This leads to the following decomposition 

h'M 1 



e -ay 2 /(2ri e ) 



h £ {y) fV 1 e-'*° 2 /(^e)o £ ( s ) s 3 £ 3fi-i ds jy^ e -as 2 /(2 Ve ) fa 

Now, let us notice that for any y £ (—1,1], 

r^ e -a S 7(2 % )0 e(s)s 3 £ 3/?-l ds 



l + e £ (y)y 6 e- 



3.3/3-1 



with D independent of e. Consequently 



< D. 



.3/3-1 



J-» ie -as 2 /(277 S )6) e ( s ) s 3 e 3/3-l ds 

1 ~^ jy^-'xs 2 /^) ds 



1 - XMe^- 1 



and there exists a constant E 1 , independent of e, such that 

sup \X e {y)\ < E. 
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Thus we can write 

K(y) .., a i-i ' 



Ml/) 



i + u £ { y y 



-m/ 2 /(2%) 



2 /(2»? £ ) ds 



and there exists a constant F, independent of e, such that 

sup \v e {y)\ < F- 
!/e(-i,i] 

Finally, for all y G (—1, 1], we have obtained 



K(y) a'M 



h e (y) g £ {y) 



< Fe 



3/3-1 



-ay 2 /(2r/ £ ) 



J\e- as2 l^)ds'' 



and the goal is now to upper-bound the last term in this equation. In this aim, we 
first consider the case where y G [— 1 + e K , 0]. In the integral, we make the change 
of variable 



V = £ ' X 



y 



2% 



with 7 = (5/3 - 2)/2, so that 7 > as soon as (3 > 2/5. We get 



-1 



where 



h{y) :-- 



^/2r] £ e^v + y 2 

Since rj £ e" / = e?/ 2 and y G [—1 + £ K ,0], with k, = (3/2, it is clear that for e small 
enough, one has: Vy 6 [— 1 + e K , 1], 



/ e (y) > £ := 



1/2 



: (1V 



'0 yJle^lH + 1 2 ' 
so that for e small enough, one has I £ (y) > 1/4. Putting all things together gives 

|A £ (y)| < ^e 3 ^- 1 x e"^ 2 = AFe* — ■+ 0, 



and the uniform convergence is proved for y G [— 1 + e K , 0]. In order to conclude 
for y G [—1 + e K , 1], it remains to notice that if y G [0, 1], one has 



-ay 2 /(2rj £ ) 



< 



1 



2a 



jy^ e - as y(2 V£ ) ds - jo^ e _ QS 2 /(2% ) ds £ ^ y 7T7]£ - 

Coming back to equation (23) gives, for all y G [0, 1] and for e small enough, 
|A £ (y)| < Fe 3 ?- 1 x eP~ l l 2 = Fe 4 ^ 2 ► 0, 

since j3 > 2/5. This concludes the case where y G [— 1 + e K , 1]. 
(2) y G (—1, —1 + s K \. let us denote y = — 1 + pe K , with < p < 1, so that our goal is 
now to upper-bound 

h' £ (-l+pe K ) 



|A e (p)| := 



<7 £ ( 



-1 



/i e (-l+pe K ) y^-l+pe*) 



that is to say 

|A £ (p)| 



e V{b E {-l+pe K ))/e e - a (-l+pe") 2 /(2 VE ) 
j-l +P e- eV(bes)/£ds ~ j-l +p£ K e _ as2/(2ve)ds 
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independently of p G (0, 1]. For any smooth function / on [—1,0], we may write 
the following Taylor expansions 

/(-l +pe K ) = /(-l) + f'(-l)pe K + f^p 2 e 2K 

and 

M ds = f(-l)pe" + ^p 2 e 2 ^ + l^Mp^ 

where 9\ and #2 belong to the interval (—1, —1 +e K ), and depend on p and e. This 
leads to 

£/■ i , k\ i 1-1- f'(~ l ) r, r K J- 1 JZ r 2n 

/(-l+pe K ) 1 +2/(-i)^ e 



Considering /(s) = e _as2// ( 2,?e \ we get 

= ap^/V" 2 " = ape 1 



/(-I) 
and 

/"(#) 2 2k _ 2( <^f_ _ a\ 2k a(l-0 2 )/(2r, e ) 

Now, since there exists q G (0, 1) such that 9 = — 1 + qe K , it is readily seen that 

e c 



e^O £^0 

and 

Thus we have the following Taylor expansion 
/(— 1 + pe K ) 1 / ap 



f:? +p£K f(s)ds Pe 

With SUPo^! \4> £ {P)\ < 00 ■ 

Now, considering /(s) = e^( 6eS )/ e , we get 

J^ve K = pe- x = ap£ 7 + &(p)e r^ 

and 

Hg) p 2 £ 2. = p 2 I ( b £ V'(b £ 0) \ 2 + gVW) \ £ 2k (V{b e e)-V(-b e ))/e 



/(-I) 

For the same reason as above, we have then 

Since /3 > 7, we have the following Taylor expansion 
/(-l+pe K ) 



-^x(l + ^ + ^(p)p e 2 ^, 



.r! 1+p£K /(a) ds Pe" 
with sup 0<p<1 |Ve(p)| < 00. Gathering the intermediate results, we get 

|A e (p)| = | Ve (p) - ^(p)|e 2 ^ K < Ge^" 4 )/ 2 , 

where G is independent of e. It turns out that the uniform convergence is ensured 
as soon as f3 > 4/9. 
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This concludes the proof of Lemma 4.10 



□ 



Step 6. The difference of the two drifts in (22) is negligible with respect to the variance 
rj e of the Brownian component as soon as 4/9 < f3 < 1/2. Theorem 4.3 in |8] ensures that 
C(Y.\T\ < T_i) and C(Z. \T\ < X!_i) are then asymptotically equivalent. This approxi- 
mation result relies on the Girsanov Theorem. The Novikov condition ensuring that the 
exponential martingale is uniformly integrable can be checked as in [5]. In particular, 

Tf ~ Tf . 

Step 7. After an obvious scaling, we have to estimate the reactive time for a repulsive 
Ornstein-Uhlenbeck process between — b e and b £ starting at 

-> 0. 



Since b £ /y/e and c e j \fe both go to infinity when e goes to zero, the estimates in the proof 
of Theorem 3.3 (see also Remark 3.5) ensure that 

Tf ~ -(-loge + log6 e + logc e + G-loga) 
e-s>0 a 

where the law of G is a standard Gumbel distribution. Putting all things together, we 
have established the following estimate. 

Proposition 4.11. Conditionally to the event {T^ < Ta], 

1 



, ->6 e ~ -(-loge + logc e + log6 £ + G-loga) 



where the law of G is a standard Gumbel distribution. 



4.4. Conclusion. The estimates of Propositions 4.1 and 4.11 are the key points of the 



proof of Theorem 1.4 



Proof of Theorem 1.4 One can write 

" B ds _ 1 
V'(s) ~ a 



B 



a 1, , 



B 



Thanks to Assumption (Q, s h-> aV'(s) 1 + s 1 is integrable on (0, B). 

1' 



Similarly, 



log h \ogB _ 1 
a a a J 



logc £ loglxl 1 
a a a 



B 



a 



a 



v , {s) +-)ds + o £ (l). 



As a conclusion, Propositions 4.1 and 4.11 ensure that, for any x £ (A,0), we have, 
conditionally to {Tg < Ta} that 



Tr 



x-±B 



b £ — vB 



1 



^ + - (- logs + log c £ + log b £ - log a + G) + i 6e ^. B 
a 

loge + log(|x|J3) + F(x) + F(B) - log a + G). 



Notice that one can let x go to A in this expression. 



□ 
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FIGURE 3. Mean time of the reactive path for the potential V given in ([6| 
as a function of loge. The 95% confidence intervals are of the size of the 
points. These results have been obtained with the algorithm described in 
[6]. The theoretical asymptotic behavior (when e goes to 0) is drawn in 
dotted line. 



Figure [3] illustrates this result for the process 

dx[ e) = -V'{xi £) ) dt + V2~edB t , 

with V(x) = x A /A - x 2 /2, X { Q £) = x = -0.89, on the set [A,B] = [-0.9,0.9]. Denoting 
7"-0.89->0.9 the length of the reactive path from -0.89 to 0.9, then Theorem |1.4| ensures 
that E[T_ .89^o.9] is equivalent to - loge + log(0.89 x 0.9) - \ log(l - 0.89 2 ) ^ log(l - 
0.9 2 ) +7 (where 7 stands for the Euler's constant) when e goes to zero. Figure [3] compares 
this theoretical result (continuous line) with the empirical means obtained thanks to the 
algorithm described in [6] for e ranging from 0.007 to 1 (circles). 

5. Other examples 

The aim of this section is to analyze the distribution of the lengths of the reactive 
paths, when the potential V has a maximum at point z* = 0, but does not satisfy the 
Assumption 1 1.1 1 More precisely, we successively consider three cases: 

(1) V behaves like — |x| around x = 0, 

(2) V is constant equal to around x = 0, 

(3) V is regular at but V"(Q) = 0. 

We will consider special potentials, for which one can derive an explicit expression for 
the asymptotic of the distribution of the lengths the reactive paths. We will see that the 



asymptotic behavior is very different from what we obtained in Theorem 1.4 



5.1. Brownian motion with drift. The easiest case to deal with is the one of the singular 
potential V(x) = — f3\x\. It corresponds to a Brownian motion with a piecewise constant 
drift, namely: 

dX^ £) = $ sgn (x ( t £) ) dt + V2edB t , 



where /3 is a positive real number and sgn(x) stands for the sign of x. In that case, 
Equation Q is a second order ordinary differential equation with constant coefficients 

Let us recall the expression of the Laplace transform of the conditionned first exit time 
on (0,6) for a Brownian motion with drift (see [3l p.309]). 
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Proposition 5.1. Choose a < x < b and and consider the process defined by 

= [it + Wt- Let us denote by H the first exit time of (a, b). Then, 

E L-sH\ w <g) = h \ _ sinh((6 - a)\fi\) sinh((x - a)y/2s + fj?) 



smh((x - a)|/i|) s inh((6 - a)y / 2s + n 2 ) 



A few remarks are in order. 



Remark 5.2. Notice that the law of H knowing that Tf, < T a does not depend on the sign 
of the drift fj,. This may seem surprising but it is consistent with the fact that going up is 
equivalent to going down after introducing the h-trans formed process, see Section 4-2 above. 

Remark 5.3. Notice that 

(24) ]imE x (e-^\W^ = b) = s j^tl^M 

V H ) M sinh((6- a)y/2s + /x 2 ) 

Remark 5.4. If \i > then H converges to Hf, the hitting time ofb as a — > —oo: 
lim E x (e- sH \W^ =b]= ecWd-^) 

which is the Laplace transform of the inverse Gaussian distribution with parameter m = 
(b — x)/fi and I = m 2 . We recall that the density of the inverse Gaussian distribution with 

parameters (to, I) is x H> \f]^x~ 3/2 exp (- '^m^c ) lx >0- 

We can use these results to study the law of the hitting of starting from x = — 5 if the 
process X& satisfies, at least when X^ 6 (—5,0): 

X ( t £) =x + \[2eB t - pt. 
From the scaling property of the Brownian motion, we can compute the Laplace transform 
F of H = inf jt > : X^ (-5,0)} conditionally to ix$ = o}, using (24): 

smh(<5/3/(2e)) y / P 2 /(2e) 2 + s/e 



FJs) 



sin 



h(5 v //3 2 /(2e) 2 + s/e) /V(2e) 



( g (l - V 1 + ^)) - ex P ("I I 1 + V 1 + W J J /., Ir.s 



For a fixed s, we thus get lim e _»o ^e( s ) = ex P ( — 7?) ' anc ^ 

4^-^))^(f>/vi> 



5s 



exp v^ 1 _ eKp ( f JT^p\ V 1+ fl- 



As a consequence, 



H > — and ■=!- > N 



e^O (3 y/e e^O \ ' ft 
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In this case, with the same reasoning as in Section 4.1 one can deduce that the length 
of the reactive path between points —6 and +6 has the deterministic limit 26 //3 when e 
tends to zero. The absence of any asymptotic randomness in the length of the reactive 
path, in contrast with Theorem 1.4 is due to the fact that in this case, we do not have 
V'(0) = 0. The next situation that we propose to deal with is the opposite one, specifically 
when V'(x) = in a neighborhood of 0, and we call it the totally flat potential. 



5.2. Totally flat potential. Let us investigate in this section the case when the potential 
V is flat around the saddle point. More precisely, let us consider the process given by 
x[ £) = V2eB t , b > and 

= inf{i>0, X[ £) i (-6,6)}. 

One has, for any s > 0, 



sinh ( y / 46 2 s/e 



Moreover, 



Notice that, for any s S 



0. W £ 



E- b (e 8H \X^ = &) = where G ( x ) 



v (-*)» ' 

l^k>0 (2fc+l)! 
In particular, 

Lemma 5.5. For any s > and b > 0, one has, conditionally to Xq = x and T5 < T_5, 
and in the limit x — > —6, 



H ie) _hU2_ 2V2 \ 

- b > b ~ e ^3 + 3^5 ) 
where E(Y) = 0, V(Y") = 1 and its Laplace transform is given by 

-sY\ VA~s Bs , , 6^ , „ V§ 



E(e ) = -, r-e where A = — — and B 

sinh(v^4s) v2 V2 



Remark 5.6. Thanks to the scaling property of the Brownian motion, this result is valid 
for any e > 0. 

In conclusion, in the case of a totally flat potential, the length of a reactive path goes to 
infinity at rate 1/e when e goes to zero. Again, this is different from the non-degenerate 
case of Theorem 1.4 where the length of a reactive path goes to infinity at a slower rate, 
namely log (1/e). 
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5.3. Degenerate concave potentials. Between the two extreme situations of Section 5.1 



(where V'(0) 7^ 0) and Section 5.2 (totally flat potential), the main result of this paper 
stated in Theorem |1.4| studies the length of a reactive path for a potential V which is non- 
degenerate at (also called quadratic case: V(0) = but V"(0) 7^ 0). In this last section, 
we briefly discuss some intermediate situations, when the second derivative of the potential 
V is equal to at the local maximum 0. Again, we will see that the asymptotic of the 



length of the reactive path is very different from the quadratic case of Theorem 1.4 To 
that end, we focus on monomial potentials: the potential V is given by 

x 2n+2 

V(x) = with n > 1. 

w 2n + 2 

We consider the diffusion process (X^) t>0 solution of 

(25) X[ £) = x + V2~eB t + / (X^) 2n+1 ds. 

Jo 

As will be explained below, in this case, the length of a reactive path goes to infinity at 

n — re 

rate e n+1 when e goes to zero. Notice that when n goes to infinity, e n+1 tends to 1/e, 



which is consistent with the scaling obtained in Section 5.2 for a totally flat potential. 



For convenience, we drop in the sequel the parameter e. Let us define 

n 1 b X 

t £ = e "+ 1 , a £ = e 2n + 2 and b £ = — , x £ = — 

a £ a £ 

and introduce the process (Xt) t>Q defined by 

X t = ^. 

a £ 

The process (JQ) t>0 is solution of the stochastic differential equation 

(26) X t = x £ + V2B t + f Xl n+l ds, 

Jo 

and we have that 



{T b <T^ b } = {f be <f_ 6e }, 



with obvious notation. On this event, T b = t £ Tb E . In Equation (26), the parameter e 



only appears in the boundary conditions as in Equation (15) for the Ornstein-Uhlenbeck 
process. Notice that, in the Ornstein-Uhlenbeck case (n = 0), t £ is equal to 1. As in 

the Ornstein-Uhlenbeck case, conditionally to the event "f^6 e < {Xt)t>o is still a 

Markov process starting from x £ and solution of 

e V{y) 

dY t = V2dB t + f £ {Y t )l {%>t] dt with f £ (y) = -V'(y) + 2 

J — b £ 

We now want to show that T bg , conditionally to ^Tb E < ^-b e |) has a limit in law when 
£ goes to zero. This will show that T b (conditionally to the event {T^ < T_j}) scales like 

n 

e n + 1 , which is the scaling announced above. 

The idea is to compare (Yt) t>0 to the solution (Zt) t>0 of the following equation 

(27) dZ t = V2dB t + f(Z t ) dt with f(z) = -V{z) + 2 

J —00 

The following lemma ensures that (Zt) t>0 goes to +00 in a finite (and integrable) time, 
even if it "starts from —00". 
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Lemma 5.7. If (Zt) t>0 is solution of Equation (27) starting from then it goes to 

+oo at a (random) finite time r e . Moreover, r e is integrable and it converges almost surely 
to an integrable random time when x goes to —oo: 

/+oo 
(p(+oo) - p(y))m(y) dy < +oo, 

where 

rx dy 



rv 

m(x) = exp ( / f(z)dz ) and p(x) = j exp ( — / f(z)dz ) dy 



m(y)' 



Proof of Lemma \5. 7| The result on the longtime behaviour of (Zf)t>o is a consequence 
of the behavior at infinity of the drift / given by Equation (27). For any x < 0, three 
successive integrations by parts lead to 

(28) ^(r-*£\<r ."«■> * < =S fi - 5^ + P " + 1)<4 - + 31 



™2n+2 I — I — ^,2n+l \ x 2n+2 1 ™4n+4 

i 

As a by-product, we get that for any x < —(2n + 1) 2 ™+ 2 , 

(29) < -x 2n+1 < f(x) < -x 2n+1 ' 2 1^ 

Let us introduce, for any n > 1, 



i _ 2n+l 

1 x 2n+2 



-oo r+oo _ g 2n+2 

C n = / e v W = / e ^+T- ds. 



oo J — oo 



For any x > 0, we have 



e 

-oo J — oo 

so that the previous computations imply that for any x > sufficiently large so that 
< C n , we have 

2n+l ^ \ ^ J2n+\ ze 



(30) < x zn+1 < f{x) < x zn+1 + 



„ _ e V{x) 
*-^n T 2n + 1 



A quick inspection of the estimates ( 29 ) and ( 30 ) indicates in particular that 



f(x) ~ |s 

| £E | — >+0O 



2n+l 



As a consequence, the process (Zf)t>o starting from i£l explodes with probability 1 at 
a (random) finite time r e and Zt — > +oo as t — > r e (see for instance [131 ch.6]. In short, 
this is a straightforward consequence of the expression of K x (T a A T^) that can be found in 
|13|. ch.6] and the fact that l/f(x) is integrable at ±oo. Indeed, for any x G (a, b), 

MTa A T h ) = - f\p(x) - p{y))m{y) dy + p( * } ~ p(a) /%(&) - p(y))m(y) 

with 

dy 



m(x) = exp / f{z) dz and p(x) 



m(y)' 

One has obviously that p(b) —> p(+oo) E (0, +oo) as 6 — )• +oo, and p(a) — > — oo as 
a — >• — oo. Thus, 

lim ~TT\ TT / (P( b ) ~ P(y)) m (y) dy = / (p(+oo) — p(y))m(y) dy G (0, +oo] . 



o-^-oo — p(a 



OG 
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Now, to show that r e is integrable (including in the limit x — > — oo), we need to prove 
that 

'■+00 

(p(+oo) — p(y))m(y) dy < +00. 

In this aim, let us first notice that for any real number y, we have 

f+00 

(p(+oo) -p{y))m{y) 



r+00 f f x \ \ ( f y 

' exp I — / f(s) ds\dx\ x exp I / f(s) ds 

/+00 / rx \ 

exp I — / f(s) ds I dx. 



Now, from the definition of /, one has clearly f(s) > s 2n+1 for any s£l. Hence for any 

y>o, 

y 2n + 2 r + OC _ x 2n + 2 

< (p(+oo) - p(y))m(y) < e~2^w I e 2 ™+ 2 dx. 



The symmetry of the potential V and an integration by parts show that for any y > 0, 

3 _a;2n + 2 / J/ „ o„_i_o 

e 2 "+ 2 dx= e v( - s > ds < 



L 



'y 

so that 



y2n+l 



< (p(+oo) -p(y))m(y) < 



r 

Since n > 1, the integrability of the function y 1— >■ (p(+oo) — p{y))m{y) when y tends to 
+00 is established. In order to conclude, we have to estimate this quantity when y goes to 
—00 as well. For this, let us first recall that 

\2n+l 

s— >— OO 



f( s ) ~ \ s \ 



so that p(+oo)m(y) is clearly integrable when y goes to —00. The estimation of the 
remaining term is slightly more involved. We rewrite it as follows 

rv 

-p{y)m(y) = -m(y) / exp(-F(x)) dx, 







where for any real number x, we define F(x) as the primitive of / with value at 

F(x) = [ X f(s)ds. 
Jo 

Notice that lim a; _ 5> _ 00 F(x) = —00. Then an integration by parts gives 

[ y < v< \\a Cn e M~F(y) f y f'(x) 
Jo exp (~ F ( x )> dx = "4 J{y) J Q J(xfi e M- F (x)) d x- 

Next, we focus on the last term of this equation, namely 

f y f'(x) 

/ TT~Y2 exp(-F(x))o?x. 

Jo f{ x r 

For this, we first deduce from the definition of / that 



V(x) / p V(x) 
f'(x) = -V"(x) + 2 rx v( . , V'(x) 



f x e v (°) ds \ v ' f x e v (*) ds i ' 
j— 00 \ J— 00 / 



From Equation (28), we know that 



P V(x) 

~ |xl 2 - +1 , 



f e y ( s ) ds x- 

J — OO 
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and more precisely that 

V U X) C V(X) „ -(2n+l) 

r x e V(s) ^ g x ->-oc \ x \ 

J — oo 1 1 



This leads to 
and 



2n 



f'(x) ~ -(2ra+l)x 



fix) , „, -(2n + l) , , . 

• ; exp(-F(x)) ~ \ 2 . 2 1 exp(-F(x)). 



From this we deduce that 

fO 



/ r^ expf— F(x))dx = — oo. 

J-oo fix) 2 



Since exp(— F(x)) = o(exp(— -F(x))) when x tends to — oo, we have 



/M 2 

rv fix) 



o fix) 



exp(— F(x)) dx = oil exp(— F(x)) dx 



and coming back to Equation (31 ) gives the following asymptotics 

" < t?( WA -exp(-F(y)) 

exp(-F(x))dx ~ —- , 

v-+-°° f{y) 



o 



so that 

rv ii 

-m(y) / exp(-F(x)) ~ —— ~ 

To sum up, we have shown that 



/+oo 
(p(+oo) - p(y))m(y) dy < +oo. 
-oo 



This ensures that 

E x (r e )= lim E x (T a /\T b ) 



a— >— oo 
i>- 



»+oo 

(p(sc) - piy))m(y) dy + I (p(+oo) - p(y))m(y) dy. 



In particular, E x (T e ) is finite for any i£i Finally, by monotone convergence theorem, 
r e has a limit almost surely when x — > — oo and 

/+oo 
(p(+oo) - p(y))m(y) dy < +oo. 
-oo 



This concludes the proof of Lemma 5.7 □ 



Thanks to Lemma 



5.7 



we see that • converges almost surely to a positive and 
integrable random variable as a — > — oo and b — > +oo. Moreover, 

/+oo 
(p(+oo)-p(j/))m(y) dy<+oo. 
-oo 

Now, notice that the drift / e that drives y is greater than /. This ensures that if Zq = Yq 
then, almost surely, Z% < It, for any t £ [0, T& e ). As a consequence, for any x G (— b E ,b), 
one has Tj^ b < T^ b . By monotone convergence, ^ converges to a random variable 
which is integrable since 

E(Tl^ bg ) < E(T£) < oo. 
To prove this result with full details, one would need to cut reactive trajectories into pieces, 
as done in Section [4] above for the quadratic case. This concludes the proof of the fact 
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that Tj) s , conditionally to < T—b s |, has a limit in law when e goes to zero, and con- 
sequently, that TJ, (conditionally to the event {T^ < Tl;,}) scales like e n+i , 

Acknowledgments. FM thanks the ASPI team of INRIA for its hospitality. 
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